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We study linear perturbations about a dust filled Bianchi Type I model with the vorticity set to 
zero. In comparison to linear perturbations about FLRW models, modes of perturbations about 
Bianchi type I models are coupled. We find that the tensor that represents the background shear 
needs to be degenerate in order for the scalar-mode perturbations to decouple from the rest of the 
flow. 
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What we currently know about our physical universe 
is based on the analysis of a model that is expanding, 
isotropic and homogeneous, and that has a cosmological 
constant. This, together with the analysis of its linear 
perturbations, undergird our picture of the universe [1]. 
The model accounts for the late time universe as is evi¬ 
dent in the analysis of large scale cosmic microwave back¬ 
ground observations [1, 2]. Although parameter determi¬ 
nation from the analysis of CMB fluctuations appears to 
confirm this picture [10], further analyses indicates that 
there are still anomalies [1]-[13]. In particular, it would 
appear that the universe may have a preferred direction 
[12]. It is not yet clear if the directional preference is 
intrinsic to the underlying model, and if so what impli¬ 
cations this would has on modern cosmology. Recently 
released Planck results [11] indicate that these anomalies 
persists. 

It was shown that some of the anomalies could be 
eliminated when WMAP first-year data was corrected 
with Bianchi IVh template [13]. The authors also find 
that the cosmological parameters required to reproduce 
the Bianchi morphology are inconsistent with what they 
called the ’’cosmic concordance” of other data sets and 
concluded that Bianchi models could not be a realistic 
physical explanation for these anomalies. What is curi¬ 
ous about these findings is that Bianchi Models appear 
ruled out because they are incompatible with isotropic in¬ 
flationary models and by the ’cosmic no-hair conjecture’ 
[14]. We need to point out that the proof of this con¬ 
jecture, a conjecture that forms the basis for the above 
analysis, assumes the weak energy condition. With the 
advent of inflationary models, the violation of the energy 
conditions opens up other possibilities. In particular, it 
is worth investigating anisotropic inflationary models [15] 
as they may address the incompatibility between Bianchi 
models and standard inflationary models. Given that the 
proof of cosmic conjecture subject to weak energy con¬ 
dition and the anomaly in the WMAP data, it may be 
too soon to discard Bianchi models. This is the motiva¬ 
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tion for this letter. We seek to answer a slightly different 
question and one that is related to the Cauchy prob¬ 
lem in General Relativity. The problem has to do with 
the preservation of the constraints equations subject to 
the propagation equations, recalling that the solutions 
to the constraints equations provide us with the initial 
conditions for the set of evolution equations. 

We use the covariant approach which differs from the 
standard metric perturbations approach in the sense that 
one does not perturb the metric and but a set of covari¬ 
ant nonlinear equations which are then linearized about 
a background of choice. The covariant equations are the 
result of splitting Einstein-Ricci-Bianchi equations, co- 
variantly, into propagation and constraint equations [18]. 
The method has been applied to studies of different kinds 
of perturbations in an almost-FLRW model [17]. In this 
letter we want to apply it to density perturbations about 
Bianchi Type I model with dust equation of state, where 
the vorticity vanishes. Previous studies of perturbations 
of anisotropic models using alternative approaches in¬ 
clude [16, 19-21, 27]. We have to address the following 
issue: immediately confronted with a major problem: 

• How does one characterize decoupled density per¬ 
turbations about Bianchi I model given that the 
scalar, vector and tensor modes are coupled for per¬ 
turbations about this background ? 

• How do we know that the conditions governing our 
characterization are preserved? 


In order to answer these questions we need to define the 
scope of the problem at hand and to identify parameters 
that play a role in the problem. 

For a general spacetime given by a well defined mani¬ 
fold and its accompanying metric ( A1 , g) , one can assume 
that there is a well-defined preferred motion of matter 
which gives rise to a unique 4-velocity. This velocity 
may be given by the average motion of matter for a fam¬ 
ily of preferred. Mathematically, the 4-velocity, which we 
denoted by (= where r is the proper time mea¬ 
sured along the fundamental world-lines. The 4-velocity 
has the property UqU® = — 1 [? ]. Physically, this velocity 
may be defined by the vanishing of the dipole of the cos¬ 
mic microwave background radiation (CMB). Only one 
4-velocity will set this dipole to zero [? ]. Using the 
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4-velocity, we can then define, and make use of the fol¬ 


lowing projection tensors. Parallel projector: 

U°-h = -u‘^Ub, (1) 

and with the properties 

U\U% = U\, U\ = 1, UabU^ = Ua, (2) 
Orthogonal projector: 

hab — Qab (d) 

which has the properties; 

h\h% = h\, h\ = 3, habU^ = 0. (4) 

We also make use of the permutation tensor Sabc, which 
is given by 

^abc — ^ ^dabc ^abc — V[abc]T ^abc'^ — 0 ; ( 5 ) 

and where rjabcd is a 4 - dimensional volume element in 
the rest-space of the fundamental observer. We denote 
the covariant time derivative of the tensor by 

= U^^eT^Kd. (6) 


This represents derivative along the observer word-line. 
The orthogonally projected covariant derivative will be 
denoted by Vq, so that the orthogonal projection a 
takes the form 

where there is total projection on all indices. If the vor- 
ticity is set to zero, V coincides with the 3-dimensional 
covariant derivatives. Using these projection tensors, 
covariant derivatives, and symmetric properties, the 4- 
velocity Ua may be split into the following irreducible 
parts: 

VaUf, = -UaUb+'^aUb=^Qhab + Crab+OJab, (7) 

where the trace 0 = VaU“ is the volume-rate of expan¬ 
sion. This rate also determines the Hubble parameter; 
H = b0. (Tab = {aUb) is the trace-free symmetric rate 
of shear tensor such that aabU^ = = 0. uiab = VjaUb] 

is the skew-symmetric vorticity tensor. 

Orthogonal projection allows the above covariant 
quantities to represent respective behaviour in the rest- 
frame of the fundamental observer. This rest frame co¬ 
incides with the orthogonal hyper-surfaces in the case 
where vorticity is absent. In general the Einstein-Ricci- 
Bianchi equations may be covariantly split into propaga¬ 
tion and constraint equations. The full non-linear set of 


equations include, the propagation equations: 

p + Qp = 0,(8) 
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dab T ^c{a^b) “f ^ab — b(10) 

Eab ^Eab CUr/^Iab ^^c{a^b) 4“ '^P^ab — b(ll) 
Hab + QHab - 3ac{aHb'j + curl{E)ab = 0(12) 


and the constraints: 



VVab = 

2 ~ 

-3Va0, 

(13) 

curl{a)ab = 

Hab: 

(14) 

V^Fab = 

^VaP + SabcfT^H^^ 

(15) 

V'TLab = 

-eabc(r\E-^. 

(16) 


We linearize these equations about a dust filled Bianchi 
I background with zero vorticity. The background quan¬ 
tities are the rate of expansion ( 0), the shear tensor 
((Tab), the energy density (p) and the electric part of 
the Weyl tensor (Eab)- Scalar-mode of perturbations are 
given by the conditions: 

VV^O,V'’0^O,V'’(TVO,i7ab = O. (17) 

The possibility of information exchange via gravitational 
waves or sound is eliminated by setting Hab = 0 and 
p = 0. Models with these conditions are called Silent 
universes [23]-[24]. The vanishing of the magnetic part of 
the Weyl tensor leads to a new constraint, curlEab = 0. 
These conditions modify the propagation and the con¬ 
straints equations as follows: 

p = -0p, (18) 

e = (19) 
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dab — ^^^ab ^c{a^b) ^abi (^9) 

Eab = -0^^ab + 3(Tc(aEb)^ - ipCTab. (21) 

and 


c\ 

= V^lTab - ^Va0 = 0, 

(22) 

C\b 

= curlffab = 0, 

(23) 

c\ 

= V''£’ab - = 0, 

(24) 

c\ 

= EabcfT^F"'' = 0, 

(25) 

^ ab 

= curlEab = 0. 

(26) 


The last constraint equation is as a result of setting 
Hab = 0. We have adopted the notation = 0 used 
in [25] where, 

= {VVab - ^Va0, curlaab - Hab, ■•} (27) 

and A= 1...5. The evolution of along is given by 

the system of equations , where do 

not have time derivatives. Detailed constraints analysis 
in the '1-1-3' formalism was given in[25] for the case of 
general nonlinear perturbations in FLRW models with 
dust equation of state, and in [26] for barotropic perfect 
fluid case for FLRW. In our analysis we find, 
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C\ + QC\ + C\-2eaVC\d = 0, (28) 

C^ab + QC\b + e^\a<^b)C\ = 0, (29) 

- ^EjC\ - ^aa^C\ + ^pC\ - ^curlC\ = Cf \ (30) 

CK + \qC\ = Ci") (31) 

C\b + \qC\^ - \c\^aEbf - + ]^pC\b = C^ab- (32) 

where the terms on the right hand side are new constraint conditions given by 

= -ea’’%Ei,'^C\d - cr\C\d) (33) 

= -SabcoJ'^a'^^^E'^d - eabcCT\ae^^E<^^- (34) 

= 3curlialE,y) - a^cS‘^\a^\e\Eb)d - (aEb)dyc - le^\a<Tb)d^eE^,. (35) 


Equation (28) requires (19), (20), (22), (23), (24), (40), (46) and (43). Equation (29) requires (20), (22), (23), (42) 
and (48). Equation (30) requires (18), (20), (21), (22), (24), (25), (40) and (46). Equation (31) requires (20), (21), 
and (25). Equation (32) requires (48), (21), (23), (25) and (26). The point about this is that the original constraints 
are said to be preserved if we can demonstrate that all the new constraint vanish. 

Using equation (20), can be written in the form given below. Equivalently, using equation (20), the electric 

(5) 

part of the Weyl tensor can be eliminated from , which leads to the form given below: 



It can be shown that the terms in equation 36 cancel 
each other out. The fully expanded form of given by 
equation 37 equivalent to equation (15) in [26], 

where it is shown there that the term vanishes when the 
shear is diagonal and degenerate. We note that our equa¬ 
tion 25 shows that one can choose a common eigen-frame 
for both shear and the electric part of the Weyl tensor, 
and one that leads to a simultaneous diagonalization of 
both. Such an eigen-frame is Eermi transported along 
the 4-velocity. This can be achieved by setting the off- 
diagonal terms of the shear and the electric part of the 
Weyl tensor set to zero. It is then possible to choose a 
tracefree-adapted irreducible frame components for aab 
and Eab defined by 
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((722 + 0-33), CT_ 


^(^22-^33), 


and where vorticity is zero. £'+/- is defined in a similar 
manner. now takes a form similar to equations 67-69 


in [26]. From our equation (25), ((t_ = 0 ^ E- = 0). 
It is then straightforward to show, using the tetrad for¬ 
malism, that cr_ = 0 (j.e.CT 22 = ( 733 ) leads to = 0. 
This means that, subject to the preceding discussions, 
all the new constraints vanish leading the preservation 
of the original constraints. This implies that the con¬ 
ditions given above for pure, or decoupled, density per¬ 
turbations about a dust filled Bianchi I model with zero 
vorticity represent consistent characterization only when 
the background shear tensor is degenerate. 

It is important to point out that this is not a proof that 
the solutions to the constraints exist only that if solu¬ 
tions exist, then they evolve consistently. The question 
of existence is one that needs to be investigated. Such 
an investigation should show how to construct a metric 
from given initial data in the covariant formalism. This 
notwithstanding, one can now study the decoupled per¬ 
turbations as will be presented in [28]. 

The following commutation relations for a scalar quan- 
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tity, / and PSTF-tensor quantities, Tab and Vab with van- ishing vorticity have been used in our analysis; 


= 2\h‘^^ah\], (38) 

eabcT^TPgV^^i = -TabS^^'^Tj’Vdp, (39) 

(Va/)' = Vaf-leVaf-aa^Vbf, (40) 

V[aV,]/ = 0, (41) 

curliT^)ab = £cd(aV^{r/)r‘'e}, (42) 

eabcT^CUrlV^’^ = T^^VaVbc - T'^^VbVac - 

(43) 

curlifTab) = fcurliT)ab + SadiaT^)V‘^f, (44) 

(VaTfcc)- = ^atbc - leVaTbc " + 2Ha^^Eae^bT^y (45) 

{V^Tab)- = ^^fab - ^0V'T,6 - a^^VeTafc + EabcH^dT^^, (46) 

(V^curlTab) = ^SabcV^^dT^^ + eabcT^d{lQ<J^^ - E-’^) - aabS^^^^aaeT^d, (47) 

(curlTab)' = curl{f)ab-^OcurlTab-<Te''Scd(a^''Tb)'^ + 3Hc{aTbf, (48) 

CUrlcUrl{T)ab = -^'^Tab + ^V(aV°r6)c + {p- ^0^)ra6 + 

STaia{Eb)^ - ^0^6)"} + aadT^‘^aab " T^’^a^a^bd + <7^‘'<Jc{aTb)d- (49) 
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